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Abstract 

Many of the mathematical models used in quasicrystal physics are based on tilings of the plane 
or space obtained by using strip projection method in a superspace of dimension four, five or six. 
We present some mathematical results which allow one to use this very elegant method in spaces of 
dimension much higher and to generate directly quasiperiodic packings of multi-shell clusters. We 
show that in the case of a two-dimensional (resp. three-dimensional) cluster we have to compute 
only determinants of order three (resp. four), independently of the dimension of the superspace 
we use. The computer program based on our mathematical results is very efficient. For example, 
we can easily generate quasiperiodic packings of three-shell icosahedral clusters (icosahedron + 
dodecahedron + icosidodecahedron) by using strip projection method in a 31-dimensional space 
(hundreds of points are obtained in a few minutes on a personal computer). 
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Quasicrystals are materials with perfect long-range order, but with no three-dimensional 
translational periodicity. The discovery of these solids in the early 1980's and the chal- 
lenge to describe their structure led to a great interest in quasiperiodic sets of points. The 
diffraction image of a quasicrystal often contains a set of sharp Bragg peaks invariant under 
a finite non-crystallographic group of symmetries G, called the symmetry group of quasicrys- 
tal (in reciprocal space). In the case of quasicrystals with no translational periodicity this 
group is the icosahedral group Y and in the case of quasicrystals periodic along one direction 
(two-dimensional quasicrystals) G is one of the dihedral groups D s (octagonal quasicrystals), 
D w (decagonal quasicrystals) and D± 2 (dodecagonal quasicrystals). 

Real structure information obtained by high resolution transmission electron microscopy 
suggests us that a quasicrystal with symmetry group G can be regarded as a quasiperiodic 
packing of interpenetrating (partially occupied) copies of a well-defined G-invariant cluster 
C. From a mathematical point of view, a G-cluster is a finite union of orbits of G in a fixed 
representation of G. A mathematical algorithm for generating quasiperiodic packings of in- 
terpenetrating copies of G-clusters was proposed by author in collaboration with Jean-Louis 
Verger-Gaugry several years ago This algorithm based on strip projection method 

has been considered difficult to use since the dimension of the involved superspace is rather 
high. The mathematical results we present in the present paper simplifies the computer 
program and allow to use strip projection method in superspaces of large dimension. 

The dihedral group D2m can be defined in terms of generators and relations as 

D 2m = ( a, b I a 2m = b 2 = {abf = e) (1) 

and the relations 

/ 7T 7T 7T 7T \ 

a(a, p) = I a cos P sin — , a sin h (3 cos — I b(a, (3) = (a, —(3) (2) 

V m mm m/ 

define an R-irreducible representation in M 2 . The orbit of D2 m generated by (a, 0) G M 2 

D 2m (a,P) = { g(a,f3) | g G D 2m } = {(«,/?), a(a, (3), a 2 (a,{3), a 2m -\a,f3)} (3) 

contains 2m points (vertices of a regular polygon with 2m sides). A two-shell Z^m-cluster 
C2 is a union of two orbits C2 = I?2m(a<i, Pi) U -D2m(a2, Z^)- 

The icosahedral group Y = 235 can be defined in terms of generators and relations as 

Y = (a, b I a 5 = b 2 = {ab) 3 = e) (4) 



FIG. 1: A three-shell icosahedral cluster is a union of three orbits of Y. 



and the rotations a, b : M 3 — >■ M 3 

a(a,/3, 7 ) = (2=Ia-£/3 + ± 7 , fa +^ + ^±7, -|a + ^ + f 7 ) 

(5) 

6(a,/3,7) = (-a,-/5,7)- 

where r = (1 + v / 5)/2, generate an irreducible representation of Y in M 3 . In the case of this 
representation there are the trivial orbit Y(0, 0, 0) = {(0, 0, 0)} of length 1, the orbits 

Y(a, cut, 0) = {g(a, ar, 0) | g E Y} where a £ (0, oo) (6) 

of length 12 (vertices of a regular icosahedron) , the orbits 

Y(a,a,a) = {g(a,a,a) \ g G Y} where a G (0, oo) (7) 

of length 20 (vertices of a regular dodecahedron), the orbits 

Y(a, 0, 0) = {g(a, 0, 0) | g G Y} where a G (0, oo) (8) 

of length 30 (vertices of an icosidodecahedron), and all the other orbits are of length 60. 
The union of orbits C 3 = Y(a,ar,0) U Y((3,(3,(3) U F(7, 0, 0)} is a three-shell icosahedral 
cluster (Fig. 1). 

In order to have explicit mathematical formulae, we start by presenting our results in a 
particular case, namely, G = D w . Let 

C 2 = L>io(o;i,/3i) U D 10 (a 2 ,p 2 ) = {ui, v 2 , v w , -v u -v 2 , -u 10 } (9) 
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where 

v\ = (v n , V 2 l), V 2 = (v 12 ,V 22 ), • v w = (v lw ,v 2 w) (10) 
be a fixed two-shell -Dio-cluster, and let 

w 1 = (v u ,v 12 ,...,v lw ) and w 2 = (v 21 ,v 22 , ...,v 210 ). (11) 

From the general theory 

as 

(a direct verification is also possible) it follows that the vectors 
W\ and w 2 from R 10 have the same norm and are orthogonal 

v 2 u + vj 2 + ... + vj 10 = v 2 21 + v 2 2 + ... + v 2 210 ^ 
(w 1 ,w 2 ) = v u v 21 +v 12 v 22 + ... +v 110 v 2W = 0. 
We identify the physical space with the two-dimensional subspace 

E 2 = { aw 1 + (3w 2 I a, (5 G R } (13) 

of the superspace M 10 and denote by E 2 the orthogonal complement 

E 2 = { x E M 10 I (x, y) = for all y G E 2 }. (14) 

The orthogonal projection onto E 2 of a vector a: G M 10 is the vector 

7T X = (x,—\ — + (x,—\— (15) 
\ K / K \ K I K 

where k = \ \wi\ \ = \ \w 2 \\, and the orthogonal projector corresponding to E 2 is 

n 1 - : R 10 — >E£ n ± x = x-nx. (16) 

If we describe E 2 by using the orthogonal basis {k~ 2 Wi, k~ 2 w 2 } then the orthogonal pro- 
jector corresponding to E 2 is 

V 2 :R W — > 1R 2 V 2 x=((x, Wl },(x,w 2 }). (17) 

The projection W 2 ,w = vr- L (Oi ) of the unit hypercube 

n 10 = [-0.5, 0.5] 10 = {(xi, x 2 , x w ) I - 0.5 < x, < 0.5 for all i G {1, 2, 10} }. (18) 

is a polyhedron (called the window of selecton ) in the 8-dimensional subspace i&f , and each 
7-dimensional face of W^io is the projection of a 7-dimensional face of Clio- The vectors 

e x = (1,0,0,0,0,0,0,0,0,0) 
e 2 = (0,1,0,0,0,0,0,0,0,0) 

(19) 

eio = (0,0,0,0,0,0,0,0,0,1) 



form the canonical basis of M 10 , and each 7-face of fii is parallel to seven of these vectors 
and orthogonal to three of them. There exist eight 7-faces of Qi orthogonal to the distinct 
vectors e^, e i2 , e ia , and the set 



X = (X!,X 2 , ...,X 10 ) 



Xi e {-0.5, 0.5} if i e {ii, i 2 , ^3} 
Xi = if % g {ii, i 2 , £3} 

contains one and only one point from each of them. There are 

10-9-8 
1-2-3 

sets of 8 parallel 7-faces of fi 10 , and we label them by using the elements of the set 
22,10 = {(M2, is) e ^ 3 I 1 < h < 8, k + 1 < H < 9, % 2 + 1 < h < 10 }. 



(20) 




(21) 



(22) 



In IR 3 the cross-product of two vectors v = (v x ,v y ,v z ) and w = (w x ,w y ,w z ) is a vector 
orthogonal to v and w, and can be obtained by expanding the formal determinant 



V X w 



i j k 



V X Vy V Z 



W X Wy W Z 



(23) 



where {i, j, k} is the canonical basis of M 3 . For any vector u = (u x ,u y ,u z ), the scalar 
product of u and v x w is 



u(v x w) — 

In a very similar way, a vector y orthogonal to nine vectors 



"My 



V X Vy V z 



W X Wy W Z 



(24) 



ui = (ttii,«i2,«i3, -,«iio) « e {1,2,3, ...,9} 
from M 10 can be obtained by expanding the formal determinant 

d e 2 e 3 ... ei 

Ull Un U13 ... Mi 10 

y = u 2 i u 22 u 2 3 ■■■ M210 
u 91 u 92 u 93 ... u 910 



(25) 



(26) 
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containing the vectors of the canonical basis in the first row. For any x G IR 10 , the scalar 
product of x and y is 



(27) 



For example, 
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is a vector orthogonal to the vectors e$, eio, w±, W2, and 



{x,y) 



Xi X 2 X3 
Vll V12 V13 
V21 V22 V23 



ei e 2 e 3 

Vll V12 V13 

V21 v 2 2 v 23 



(28) 



(29) 



for any x G Mr. The vector y belongs to E% , and since — 7r e$ is a linear combination of 
wi and u»2, it is also orthogonal to n^-e^, vr^es, n ± eio. Therefore, y is orthogonal to the 
7-faces of W2,io labelled by (1, 2, 3). Similar results can be obtained for any (ii, i 2 , i 3 ) G X 2i io- 
Consider the strip corresponding to W 2i i (Fig. 2) 



'2,10 



{x G 



>io 



n ± x G W 2 , w } 



(30) 
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FIG. 2: The window and the corresponding strip (case of a one-dimensional physical 

space E\ embedded into a three-dimensional superspace). 



and define 



«2«3 



max 

a 3 - £{0.5, 0.5} 



OL\ CH 2 O-Z 
Vu x Vu 2 V lia 
V2h V 2i2 V 2i3 



for each (ii, i 2 , ^3) G 2^2. 
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(31) 



A point x = (xi, x 2 , Xio) G M 10 belongs to the strip £2,10 if an d only if 



Vlh Vli 2 v li3 

V2h V 2i 2 V 2i 3 



< di 



Hll2%a for any (i h i 2 , z 3 ) G X 2 ,io- 
The set defined in terms of the strip projection method J, | 



(32) 



0.2 = V 2 {S 2;10 n z 10 ) = {v 2 x\ x g s a ,io n z lu } 



10 



(33) 



can be regarded (Fig. 3) as a quasiperiodic packing of translated (partially occupied) copies 
of C 2 . Since 



V 2 ei = {(e h wt), (e h w 2 )) = (v u ,v 2i ) = Vi for any i G {1, 2, 10} (34) 



we get 



V 2 ({ x ± ei, x ± e 2 , x ± ei } n S 2 ,io) 

C {V 2 x ± ui, P 2 x ± v 2 P 2 x ± ^10} = T^z + C 2 



(35) 



that is, the neighbours of any point V 2 x G Q 2 belong to the translated copy V 2 x + C 2 of C 2 . 



FIG. 3: The two-shell decagonal cluster C 2 and a fragment of Q 2 . 



Our algorithm works for any finite group G, any IR-irreducible representation of G in a 
space M, n and any G-cluster symmetric with respect to the origin 

C n = {v lt v 2 , v k , -v u -v 2 , -v k } C R n . (36) 

We identify the 'physical space' with the n-dimensional subspace 

E n = {a 1 w 1 + a 2 w 2 + ... + a n w n \ a ± , a 2 , a n E R } (37) 

of the superspace M. k spanned by the orthogonal vectors with the same norm 

Wi = (vn, v i2 , v ik ) i E {1,2, ...,n} (38) 

defined by using the coordinates % of the vectors Vj = (vij,v 2 j, ...,v n j). 

The window W n ^ = ^(Qk) is a polyhedron in the (k — n)-subspace E^. Each (k — n — 1)- 
face of W n ,k is the projection of a (k — n — l)-face of the unit hypercube Q, k = [—0.5, 0.5] fc . 
Each (k — n — l)-face of Q k is parallel to k — n — 1 vectors of the canonical basis {ei, e 2 , e k } 
and orthogonal to n + 1 of them. For each n + 1 distinct vectors e^, e i2 , e in+1 the number 
of (k — n — l)-faces of fl k orthogonal to them is 2 n+1 . There are k\/[(n + 1)! (k — n — 1)!] 
sets of 2 n+1 parallel (k — n — l)-faces of fl k . In the case n = 2 these sets are labelled by 

X 2 , fc = {(ii, 12,13) G Z 3 I 1 < ii < /c - 2, ix + 1 < i 2 < k - 1, i 2 + 1 < i 3 < k } (39) 

and in the case n = 3 are labelled by 



1z,k = { (h, i2,i3,U) e Z 4 



1 < ii < A; — 3, ii + 1 < i 2 < k - 2, 
i 2 + 1 < 13 < k — 1, i 3 + 1 < 14 < k 
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(40) 



FIG. 4: Projected positions of points of C 3 and of Q 3 down a fivefold axis in the case of an 
icosahedral three-shell cluster. 



In the case n = 3, a point x = (x 1 ,x 2 , x k ) € M fc belongs to the strip «S 3ifc if and only if 



Vlil ^1*2 ^3 ^4 
f2ii ^2i 2 ^2j 3 V 2i 4 
V3h V3i 2 V 3 i 3 V 3ii 



< d 



1112«3*4 



for each (i u i 2 , i 3 u) e 23,* (41) 



where 



Cti tt2 CC3 «4 

flii ^li 3 ^1*4 

V 2il V 2i2 V 2is V 2ii 
V3h V 3 i 2 V 3i3 V 3i4 

The pattern defined in terms of the strip projection method 



d-iii 2 i 3 i4 max 

ay €{0.5,0.5} 



Qn = V n (S n , k n Z fc ) = {V n X \ X E S n , k fl Z k } 



where 



V n : R k — > R n P n a; = (fowi), (x,iu 2 ), <x, «;„)), 
can be regarded as a quasiperiodic packing of copies of the starting cluster C n . 



(42) 



(43) 



(44) 
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Projected positions down a fivefold axis of points of a fragment of the pattern Q3 obtained 
by starting from the icosahedral three-shell cluster 

c ^ Y faW vfe' °) u Y (tt Ts' 7i) u m °' 0) (45) 

are presented in Fig. 4. The dimension k of the superspace in this case is 31. 

The description of the atomic structure of quasicrystals is a very difficult problem. Elser 
& Henley and Audier & Guyot jl(| have obtained models for icosahedral quasicrystals 
by decorating the Ammann rhombohedra occuring in a tiling of the 3D space defined by 
projection. In his quasi- unit cell picture Steinhardt ^| has shown that the atomic structure 
can be described entirely by using a single repeating cluster which overlaps (shares atoms 
with) neighbour clusters. The model is determined by the overlap rules and the atom deco- 
ration of the unit cell. Some important models have been obtained by Yamamoto & Hiraga 

, Gratias Puyraimond and Quiquandon [15[ by using the section 
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method in a six-dimensional superspace decorated with several polyhedra (acceptance do- 
mains). Janot and de Boissieu |l6( have shown that a model of icosahedral quasicrystal can 
be generated recursively by starting from a pseudo-Mackay cluster and using some inflation 
rules. In the case of all these models one has to add or shift some points in order to fill 
the gaps between the clusters, and one has to elliminate some points from interpenetrating 
clusters if they become too close. 

The strip projection method has been used (to our knowlege) mainly for generating tilings 
of plane or space, and only in superspaces of dimension four, five or six. Algorithms and de- 
tails concerning computer programs for generating quasiperiodic tilings have been presented 
by Conway and Knowles |7[, Vogg and Ryder [l7|,Lord, Ramakrishnan and Ranganathan 
18| . We present an algorithm for generating quasiperiodic packings of multi-shell clusters 
based on strip projection method used in a superspace of large dimension. The quasiperi- 
odic sets obtained in this way have the remarkable mathematical properties of the patterns 
obtained by projection and the desired local structure. Each point of the set is the centre of 
a more or less occupied copy of the starting cluster, and the clusters corresponding to neigh- 
bouring points share several points. Computer programs in FORTRAN 90 are available via 
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